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THE COSMIC GAMMA-HAY SPECTRUM FROM SECONDARY 
PARTICLE PRODUCTION IN COSMIC-RAY INTERACTIONS 


ABSTRACT 

The cosmic gamma-ray spectrum below 1 GeV arising from cosmic- 
ray p-p interactions is calculated. Its properties are determined by 
the properties of secondary neutral pion production occurring at ac- 
celerator energies. An "isobar-plus-fireball” model is chosen for nu- 
merical calculations in which the two dominant modes of neutral pion 
production at accelerator energies are the production of the A (1.23 8) 
isobar and one fireball. Each of these modes of production then yields 
a component to the total gamma-ray spectrum. The affect of a-p and 
p-a interactions on the cosmic gamma-ray spectrum is also calculated. 
The components are summed and the results are given in terms of both 
differential and integral gamma-ray energy spectra. 
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THE COSMIC GAMMA-RAY SPECTRUM FROM SECONDARY 
PARTICLE PRODUCTION IN COSMIC-RAY INTERACTIONS 


I. INTRODUCTION 

In previous work, the general formulas were derived for the calculation of 
the gamma-ray spectrum due to secondary particles produced in cosmic-ray 
interactions with atoms of the interstellar or intergalactic gas (Stocker, 1966). 

A calculation was then made of the gamma-ray spectrum from cosmic-ray inter- 
actions for gamma-ray energies greater than 1 GeV (Stecker, Tsuruta, and 
Fazio, 1968), It was pointed out that the gamma-rays arising from the decay of 
heavy baryon secondaries (hyperons and nucleon isobars) determine the charac- 
teristics of the gamma-ray spectrum at energies greater than 6 TeV. 

This paper deals with the detailed calculation of the cosmic gamma-ray 
spectrum from secondary particle production with particular emphasis on the 
energy region below 1 GeV. This region is of particular interest at px-esent be- 
cause the recent detection of cosmic gamma-rays in the region has indicated the 
possibility of spectral measurements in the near future (Clark, Garmire and 
Kraushaar, 1968) and because spectral measurements in this region will allow 
us to determine whether secondary particle production or another process is of 
prime importance in the production of galactic gamma-rays (Stecker, 1969). 

Several authors (Pollack and Fazio, 1963, Hayakawa, Okuda, Tanaka and 
Yamamoto, 1964, Ginzburg and Syrovatsldi, 1964) have addressed themselves to 
this problem in the past. However, in this paper we present a more detailed 
calculation based on the results of recent accelerator experiments. 

The general features of gamma-ray spectra resulting from neutral pion- 
decay should be kept in mind whenever a discussion of a particular result is de- 
sired. The general results will prove invaluable as a guide in discussing the 
results of the particular model used in the numerical calculations. 

These features may be summarized as follows (Stecker, 1966): 

A. If we plot the general gamma-ray spectrum from neutral pion decay 
on a logarithmic scale in Ey, the resultant curve will be symmetrical about 
the line in E y = in (m^/2). 
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B. The spectrum will have its maximum value at the point E y ** m n /2. 

C. Higher energy pions produce a broader spectrum than their low-energy 
counterparts. They produce the gamma rays at the high- and low-energy ex- 
tremes of the spectrum, 

A specific model will be ohosen that will be used for numerical calculation 
of the gamma-ray spectrum below 1 GeV, This model is of necessity an over- 
simplification of the true physics of pion production, but it can be made to fit 
the present accelerator data to an extent that is adequate for our purposes. The 
model is of the "isobar-plus -fireball" type. The terms "isobar" and "fireball” 
are here loosely defined to characterize the main features of the model. The 
model assumes that at accelerator energies there exist two dominant modes of 
pion production. Thus we may speak of two pion components; there is the isobar 
component, which dominates pion production in collisions where the cosmic-ray 
protons have energies in the low GeV range, and there is the fireball component, 
which supplies the main bulk of pions for p~p collisions of initial energy greater 
than 5 GeV. 


II. THE ISOBAR-PLUS-FIREBALL MODEL 


In this model, we will assume that in collisions of accelerator energies all 
pions are produced in either of two ways. 


A. Via intermediate production and decay of the A (1.238) nonstrange isobax'. 
This isobar is assumed to carry momentum directly forward or backward in the 
center~of~momentum system (c.m.s.) of the collision with an equal probability. 
The assumption made here is very similar to the one made in Steeker, et. al. 
(1968) for the other isobars and hyperons discussed. 


B. From the energy remaining available in the c.m.s. of the collision, a 
thermal pion "gas" is created where the pions are given an energy distribution 
which is very similar to a Maxwell-Boltzmann type distribution (Morrison, 1963). 
This phenomenon is sometimes referred to as a fireball in cosmic-ray 
research and will be referred to as a fireball here. Its conception dates back to 
some of Fermi’s early research into the problem of pion production in cosmic-ray 
collisions (Fermi, 1951). Our low-energy fireball is assumed to be roughly at 
rest with respect to the c.m.s. of the collision with its pions directed isotropically 
in this system. As an illustration, Figure 1 shows schematically the situation 
resulting from a TeV (10 3 GeV) collision as pictured in previous work (Steeker, 
el. al., 1568). The large difference in the behavior of the mesons and baryons 
emanating from the collision leads to the results obtained for the high-energy 
gamma-ray spectrum described in that paper. 
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The situation for low-energy collisions is not that simple, as is shown in 
Figure 2. As is indicated, it is much harder to distinguish between the two modes 
of pion production because of the similarity in the momenta involved. There is 
no longer the large disparity in c.m.s. momenta of the previous situation. 

Indeed, we may even conceptually fuse the two components and envision all 
the created pions together as having an elongated c.m.s. momentum distribution 
whose anisotropy is a function of collision energy. A simple model of this type 
has been used by Hayakawa, Okuda, Tanaka, and Yamamoto (1964) in their calcu- 
lations of the cosmic gamma-ray spectrum. However, the two-component model 
presented here utilizes the more detailed results of more recent accelerator ex- 
periments. Among the most striking of these results is the complete dominance 
of the A (1,23 8) production mode in collisions where E < 3 GeV (Muirhead, 1965 

р, 666), Indeed, even for collisions below the threshold for production of the mass 
peak at 1.238 GeV, there is evidence that this resonance (whose mass distribution 
extends below 1.238 GeV) affects the pion energy distribution (Focardi, Saporetti, 
Bertolini, Grigoletto, Peruzzo, Santangelo, and Gialanella, 1965). We will assume 
here that below 3,16 GeV, all pion production occurs through the intermediate 

A (1.238) resonance. This hypothesis is in agreement with the observation that 
at about 0,6-GeV proton kinetic energy, the neutral pion distribution is isotropic, 
since in this energy range the isobar enhancement would be formed at rest in the 

с. m.s. The model is also in agreement with the observation of an increasing 
asymmetry in pion momentum toward the forward and backward directions as 
E p increases, since at higher energies the isobars tend to carry more end more 
momentum in the forward and backward directions. 


III. CALCULATION OF THE GAMMA-HAY SPECTRUM FROM 

THE ISOBAR PION COMPONENT 

We are now ready to specify the parameters needed for the numerical calcu- 
lation of the gamma-ray spectrum from the isobar pion component, hereafter 
referred to as the i-process component. For our model we need to specify 

A. the energy distribution of the isobar in the c.m.s, 

B. the angular distribution of the isobar in the c.m.s. 

C. the energy distribution of the decay pions in the BARS (baryon-isobar 
rest system) . 

D. the angular distribution of the decay pions in the BARS. 
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E. the cross section for production of the A (1.238) irobar as a function of 
initial proton energy. 

In our model (Stecker, et. al., 1968) the isobars are assumed to decay iso- 
tropically in the baryon rest system (BARS). The energy distribution of the 
decay pions is obtained from the Breit-Wigner form of the isobar mass distribu- 
tion. The isobar mass distribution is assumed to peak at Mj » 1.238 GeV and have 
a width of 0.1 GeV. The Breit-Wigner distribution is given by 


3(M*) = i 

v (m* - Mj) 2 + r 2 


( 1 ) 


We denote the pion energy in the BARS by the symbol /a. In terms of isobar 
mass, jj. is then given by 


M* 2 + M 2 - M 2 

fi~ 

2M* 


( 2 ) 


The normalized distribution function for fi is related to the normalized 
Breit-Wigner mass distribution by the transformation 


V) 


3(M*) 

d^TdivP 


(3) 


which then yields 


30 ) =~J (m + 4?" + (M 2 - M^) 2 - Mj) 2 + r 2 ] 


1 - 




r l 


m + 4^ 2 + m 2 - 


(4) 


4 



Equation (4) was solved on the CDC 6400 computer as a check on on inter- 
mediate step in the complete numerical solution. The results are graphed in 
Figure 3 along with the BARS gamr a-ray distribution resulting from this pion 
distribution and given by 



pm 

"By + / 4Ey 


// jl 2 - M* 


(5) 


In the actual calculation of the gamma~ray spectrum, we must take into 
account the kinematic limits on the allowed pion energy. There is only a limited 
amount of energy available for pion production. In order to specify a maximum 
BARS pion energy, ^ max , we note that in the extreme case where all the c.m.s. 
energy of the collision Is transformed into the mass of the resonance we "ind 


and thus we can define 


M 


max 


(2 y 0 - *) M p , 


IA 


M * 2 + M 2 

max tt 


K 


max 


2M. 


max 


( 6 ) 


( 7 ) 


We must then renormalize the distribution function over /x so that its integral 
over the allowed energy region is equal to unity. We thus define a renormalizing 
weighting factor 


w E 

rn 



“1 


dfJL S(/x) , 


( 8 ) 


which multiplies the distribution 5(^). 
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Wo specify In our model, as before, that the isobars carry momentum either 
directly forward or directly backward in the c.m.s. (Morrison, 1983). We further 
assume that the low-energy i-process produces pions through the two-stage decay 


p + p ** p + A + (1,238) 

L P + 77° 

1 1 hi >ii 2 y ( 


(9) 


Let E G be the total energy of the collision in the c.m.s., i.e. 

E c = 2M p 7 c . (10) 

Furthermore, let us designate the energy of the isobar and proton in the c.m.s. 
system of the collision by E ci and E cp , respectively. Then 



E 2 + M* 2 - M 2 





( 11 ) 


The energy of the forwardly and backwardly produced isobars in the laboratory 
system are then 


Ef =r=(E ci +/3 C P c5 ), (12) 

and 

E b " (E ei “ P c i)> < 13 > 

with 

Pci = /E ci - M* 2 . (14) 


We may therefore write 




-E f ) + 8(E. ~E b )]. 


(15) 
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Using the results of previous calculations (Stecker, et aL, 1968), it now follows 
that the gamma-ray spectrum from the A (1,238) i-process is given by the relation 


W) 






dE„ I(EJ w rn (E„) | dE. 


S 


M* 

/E? - M ’ 2 


^(E.; E ) Rj 


^(Ey.E;) 

A(e 7 ,e 1 ) 



l 


ft (E n ) 
r raax v P' 


E y + M*./ 4E y 


dfj. 


5 (^) 

/fj? - M* 


(16) 


where w rn is defined by equation (8), /jl ax is defined by equation (7), E^, is the 
gamma-ray energy in the BAHS, = 2?3*is the branching ratio for the neutral 
pion decay mode of the A(1.238) isobar, 


and 


U(E y , E.) 


E 


7i(l ~ PiY 


L(EL, E.) = 


E 


r 


7i (1 + P L ) 


(17) 


(18) 


All the integrations in equation (16), except for the dE p integration, were 
carried out numerically by Gaussian quadrature on a CDC 6400 computer. The 
integration over dE p was carried out in the following manner: 

The product R. cr, (E p ) (in millibarns) was approximated by the discontinuous 
function. This function was taken to be a three-part power law of the form 

\ 

\ 
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0 

for 

E,<1.48 GeV 

10" 1 E 0 ’ 7 

P 

for 

1.4S < E 0 < 1.88 GeV 

r 

6.9 

for 

1.88 < E p < 3.16 GeV 

69 E~ 2 

p 

for 

3.16 < E p < 10 GeV, 



( 19 ) 


with cr. in mb.* 

1 

The differential cosmic-ray proton intensity is shown in Figure 4. 


IV. CALCULATION OF THE GAMMA-HAY SPECTRUM FROM THE 
FIREBALL COMPONENT 

We now calculate the gamma-ray spectrum produced by pions created in the 
c.m.s. of the collision. We will refer to these pions as the fireball component. 
(The use of the term here is restricted to our previous definition.) For the 
purpose of discussion, we define the f-process (as opposed to the i-process) as 
the following production and decay scheme: 


P + P - P 



(20) 


where £ ± and £ 0 are the average numbers of charged and neutral pions, re- 
spectively (a function of initial proton energy), and F is the so-called fireball. 

By no means should F be regarded as a real particle; it is better to picture it as 
a bosonic nonparticle whose only properties are 1) a momentary possession of 


*Such a cross section is an approximation which assumes complete dominance of the i-process 
for neutral pion production in collisions involving kinetic erorgies below 1 GeV but which main- 
tains the observed E~^ energy dependence of the A(l,238) cross section in various channels for 
kinetic energies above 2 GeV (Muirhead, 1965, W. E. Ellis, D, J. Miller, T. W. Morris, R. S. 
Panvini and A. M. Thorndike, 1968, Phys. Rev. Letters 2J, 697, P, L. Connolley, W. E. Ellis, 

P. V. C. Hongh, D. J. Miller, T. W. Morris, C. Onannes, R. S. Panvini and A. M. Thorndike, 1967, 
BNL internal report 11980. 
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the mass -energy of the created plons, and 2) a momentary occupation of the 
interaction volume centered at the center of momentum of the collision and 
therefore moving at p with respect to the laboratory system. 

We assume that the fireball decays isotropically with a c.m.s. energy dis- 
tribution which is roughly independent of energy. The justification for this 
assumption comes both from accelerator studies (Morrison, 1963*, Fidecaro, 
Finoochiaro, Gatti, Giacomelli, Middelkoop, and Yamagata, 1962) and from 
studies of interactions of cosmic-rays in the atmosphere (Fowler and Perkins, 
1964). Accelerator studies of*p~P interactions for energies up to 30 GeV give a 
charged-pion transverse-momentum distribution that fits an empirical law of 
the form 


{ rr(Pl)* Pl e'^^n 


(21) 


(Morrison, 1963). Studies at 300 GeV (Guseva, Dobrotin, Zelevinskaya, Kotelnikov, 
Lebedev, and Slavatinsky, 1961) fit this law very well. (There is only a change in 
normalization involved, indicating an increased multiplicity.) It has also been 
observed that is roughly constant in p-p collisions involving energies of 

up to 10 4 GeV (Fowler and Perkins, 1964). If neutral pions also obey this law, 
we would expect that the resultant gamma rays would have a transverse- momentum 
distribution of the form 


£y(Pj.) a e p ^- /<p L > t 

with 

< P l>r = I < p i\‘ 


( 22 ) 


(23) 


Data confirm this exponential distribution for > 0.1 GeV in p-p collisions 
of energies up to 10 4 GeV (Fowler and Perkins, 1964). Accelerator experiments 
(Fidecaro, et. al., 1962) mesuring the gamma-ray spectrum from p-p interactions 
at 23.1 GeV have substantiated indications of an exponential distribution of neutral 
pions in the c.m.s. The data are well represented by an energy distribution 
function of the form 



( 24 ) 
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Thus,, over a wide range of energies 


UP.!.) a Pi 


e -2p i /<p i > 7r 


(25) 


The quantity is a relativistic invariant that is proportional to the c.m.s. 
momentum of the pion, and, therefore (relativistic ally) to its c.m.s. energy. 
Letting £ stand for c.m.s. pion energy, we may therefore use the approximation 


3(f) a f e” 2 ^ /< ^ > , 



Since <(f)> ~ 0.5, independent of energy for collisions involving energies up 
to 10 4 GeV, we find that 


0(0 = ~ € e - *, 


(27) 


with the normalization constant. Z defined by 


r^mnx( E p) 

Z= (28) 


so that 



df 3(f) = 1 . 


(29) 


Evidence for a two-component pion production model with the f-process 
component producing pions isotropically in the collision c.m.s. has been pre- 
sented by Dekkers, Geibel, Mermod, Weber, Willitts, and Winter (1964). Their 
work indicates that the i-component yields a very large ratio of positive-to- 
negative pions. Thus, the negative pions arise principally through the f-process. 
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It is then found that at accelerator energies, most negative pions are produced 
isotropically in the c.m.s. (the remaining ones being consistent with a small 
contribution from N* -* N v n decays). 

We have thus assumed that the i-process plus the 1-process provide all the 
pions produced, i.e., 


Mr 2 U(E p )a(E p )] TotQl -o,(E p ). (30) 

The function ( l (E ) a(E )L t is shown in Figure 5. 

A gamma ray of energy e£ c ) in the c.m.s. may have any laboratory energy 
E y within the range 


T c E<, c > (1 - /3 o ) < E y < y c E<, c > (1 + /3.) . 


(31) 


Under our assumption of uniform c.m.s. pion emission for the f-process, there 
is a constant production probability, 


f(E y ; E< c >) = 


1 


2y c P c E^) 


(32) 


within this energy range. It follows from equation (31) that the integral over the 
source function with respect to E^ c) is bounded by the upper and lower limits 
U F and L p , respectively, where 



E 


y 


y.o-/ 3 .) 


and 


MV = 


E. 


y 


'M 1 * V’ 


(33) 


(34) 


since only gamma rays within this c.m.s. energy range will produce gamma rays 
of laboratory energy E y . 
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Therefore, the complete formula for the gamma-ray spectrum arising from 
the f-process is 


I 


f 




^(y^y i<v 

y c (E p ) p Q (E p ) 


/“°F< E p> dE< c > r^ma 

4<v ~w> *4° 




df. 


m 


; E< / C >+m^/4E^> /eZrf 


(35) 


The resulting f-process spectrum was obtained numerically, using the same 
methods previously described for the evaluation of the i-process spectrum. The 
results of the numerical integration are shown in Figure 6. Figure 6 also shows 
the total spectrum produced by summation of the i- and f-components . The integral 
spectrum is given in Figure 7. 


V. COSMIC GAMMA-HAYS PRODUCED IN COSMIC-RAY p~a AND 
a-p INTERACTIONS 

Up to this point, we have restricted ourselves to a discussion of cosmic 
gamma-ray production in cosmic-ray p-p interactions. The details of particle 
production in p-p interactions have been well studied in recent accelerator ex- 
periments. Such is not the case for interactions involving protons and helium 
nuclei. However, these interactions may be expected to provide a substantial 
contribution to the total gamma-ray production of roughly the same magnitude 
as the contribution from p-p interactions and therefore these interactions 
cannot be neglected. 

Figure 5 shows the total cross section times multiplicity for neutral pion 
production in proton-helium interactions as a function of kinetic energy per nucleon. 
Below 1 GeV/nucleon kinetic energy, these cross sections have been determined 
experimentally (Prokoshkin and Tiapkin (1957a, b), Batson, Culwick, Klepp, and 
Riddiford (1959), and Kozodaev, Kalyukin, Sulyaev, Filippov, and Shcherbakov 
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(I960)). Below 1 GeV/nuclecn kinetic energy, it has been found experimentally 
that neutral pion production precedes predominantly through the A (1.23 8) pro- 
duction channel, this manifested both by the angular distribution of the pions and 
the 2:1 ratio of neutral to positive pion production in p-n interactions (Prokoshkin 
and Tiapkin, 1957a). We have therefore made the assumptions that the same 
isobar-plus-fireball model used for p-p interactions is also valid for p-a and 
a-p interactions, that the i-process dominates below a total energy of 3.16 GeV/ 
nucleon and that the energy dependences of the cross sections for the i-process 
and the f-process are the same for nucleon- nucleon collisions in all of these 
interactions (charge independence). For interactions involving kinetic energies 
greater than 1 GeV/nucleon, the cross section data was extrapolated by assuming 
that the energy dependence of the pion multiplicity has the same form for a-p 
and p-a interactions as for p-p interactions. We assume also that the cross 
section for A(1.238) production decreases with energy as E" 2 as is the case for 
p-p interactions. The assumption is also made that the pions are produced pri- 
marily in interactions which simulate free nucleon- nucleon interactions (Batson, 
et. al. (1959)) but with shadowing corrections as given by Kozodaev, et. al. (1960). 

The energy spectrum of cosmic-ray alpha particles at solar minimum is 
shown in Figure 4. These particles produce gamma-rays by interacting with the 
interstellar gas (primarily a-p interactions). To this contribution, we must add 
the contribution of cosmic-ray protons interacting with helium nuclei in the 
interstellar gas (p-a interactions) . Interstellar helium is estimated to make up 
approximately 36% of the nucleons in the interstellar medium so that the ratio 
of helium nucleons to hydrogen nucleons is 0.57 (Allen, 1963). 

The differential energy spectrum of cosmic gamma-rays from a-p and p-a 
interactions was calculated based on the intensities and cross sections given in 
Figures 4 and 5 respectively and on the assumptions given above. The results 
of this calculation, when added to the results of the p-p calculations, yield the 
total spectrum given in Figure 6. The total integral production spectrum is 
shown in Figure 7. 
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FIGURE CAPTIONS 


Figure 1 - A schematic representation of the distribution of secondary pions 

resulting from a nucleon- nucleon interaction in the TeV energy range. 

Figure 2 - A schematic representation of the distribution of secondary pions 

resulting from a .nucleon-nucleon interaction in the GeV energy range. 

Figure 3 - The energy distribution function of pions in the BARS resulting from 
the decay of the A (1 .238) isobar and the energy distribution function 
of gamma-rays in the BARS resulting from the decay of these pions. 

Figure 4 — The differential cosmic-ray proton and alpha particle energy spectra. 

Figure 5 - Cross section times multiplicity for neutral pion production. Data 
for p-p interactions as summarized by Muirhead (1965). Data for 
a-p interactions from Prokoshkin and Tiapkin (1957a, b), Batson, 
Culwick, Klepp, and Riddiford (1959), and Kozodaev, Kulyukin, 

Sulyaev, Filippov, and Shcherbakov (1960). 

Figure 6 — The calculated differential production spectrum of gamma-rays 
produced in cosmic-ray interactions based on the isobar-plus- 
fireball model described in the text. 

Figure 7 — Calculated integral production spectrum based on the results of 
Figure 6. 
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Figure 1— A schematic representation of the distribution of secondary pions resulting from a 

nucleon-nucleon interaction in the TeV energy range. 
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Figure 2-A schematic representation of the distribution of secondary pions resulting 
from c nucleon-nuclein interaction in the GeV energy range. 


18 



0 0.200 0.400 0.600 0.000 1.000 1.200 


E(GeV) 

Figure 3-The energy distribution function of pions in the BARS resulting from the 
decay of the A (1.238) isobar and the energy distribution function of gamma-rays 
in the BARS resulting from the decay of these pions. 
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Figure 5~Cross section times multiplicity for neutral pion production. Data for p-p inter- 
actions as summarized byMuirhead (1965). Data for a-p interactions from Prokoshkin and 
Tiapkin (1957a, b), Batson, Culwick, Klepp, and Riddiford (1959), and Kozodaev, Kulyukin, 
Sulyaev, Filippov, and Shcherbakov (1960). 
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